Observations and models of lateral hydrothermal circulation on a young ridge flank: Numerical evaluation of thermal and chemical constraints

[1] We used a two-dimensional coupled heat and fluid flow model to investigate large-scale, lateral heat and fluid transport on the eastern flank of the Juan de Fuca Ridge. Cool seawater in the natural system is inferred to enter basement where it is exposed close to the spreading center and flow laterally beneath thick sediments, causing seafloor heat flow to be depressed relative to that input at the base of the plate. The flow rate, and thus the properties of permeable basement (the flow layer), controls the efficiency of lateral heat transport, as quantified through numerical modeling. We simulated forced flow in this layer by pumping water through at a fixed rate and quantified relations between flow rate, thickness of the permeable basement, and the extent of suppression of seafloor heat flow. Free flow simulations, in which fluid flow was not forced, match heat flow constraints if nonhydrostatic initial conditions are used and flow layer permeabilities are set to the high end of observed values (10 11 to 10 9 m). To match seafloor heat flow observations, the models required lateral specific discharge of 1.2 to 40 m/yr for flow layer thicknesses of 600 to 100 m, respectively. The models also replicate differences in fluid pressures in basement, and the local distribution of pressures above and below hydrostatic. Estimated lateral flow rates are 10 to 1000 greater than estimates based on radiocarbon ages of basement pore waters. Estimated lateral flow rates based on thermal and chemical constraints can be reconciled if diffusion from discrete flow zones into less permeable stagnant zones in the crust is considered.